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1. Introduction 

Clusters are a distinct state of matter. Not only their structure, electronic, magnetic 
and optical properties change with size, but also new properties emerge as a result of 
quantization effects unknown in the atom and the bulk solid [1]. Well known examples 
for unique cluster properties are magic sizes in geometric and electronic structure, 
superparamagnetism, plasmon resonances or size dependent reactivity [2l[3lll]. Though 
cluster properties may be studied in their purest form for mass selected clusters in 
the free beam, a use of clusters requires a suitable support [5]. Ideally, a system of 
supported clusters is a regular array with equal sized, equally spaced clusters, each in 
an identical environment and of macroscopic array extension. The cluster bonding to 
the substrate should be strong enough to warrant their stability at the temperature of 
operation and weak enough not to destroy them as own entities. The substrate should 
be inert in a sense that it does not deteriorate under the condition of use but active 
to bestow new functionalities to the clusters without destroying them. Such arrays 
or cluster superlattices would allow one to address single clusters (as might e.g. be 
of interest for magnetic applications), to obtain a large amplitude of response resulting 
from the additive superposition of all single clusters and being characteristic for a cluster 
of well defined size in a specific environment (as e.g. needed in catalysis), or to obtain 
a collective coherent response characteristic of a large interacting ensemble of clusters 
(as e.g. in optics). 

This vision triggered attempts to make use of self organization and/or templates for 
the creation of cluster superlattices and a few examples are given here (compare also [6]). 
By making use of the regular arrangement of steps on vicinal Au(lll) and its herringbone 
reconstruction crossing the steps at right angles a proper Co cluster superlattice could be 
realized [7] , enabling the measurement of its temperature dependent magnetic properties 
including the distribution of cluster magnetic anisotropy energies [H] • Also Fe could be 
grown regularly [9], but the method is limited by its low growth temperature and strong 
cluster substrate interaction giving rise to alloying. Alumina double layers on Ni3Al(lll) 
exhibit a phase with a regular superlattice of sites for nucleation upon metal deposition 
[To] , which are holes in the oxide [11]. Cluster superlattices were realized for a variety of 
metals at room temperature [TOl [T2l [T3] and a high perfection was achieved for Pd [TO] . 
Though such cluster lattices appear to be suited as model catalysts, the preparation of 
the alumina layer is the result of a subtle procedure, comprises defects and at least two 
phases of alumina with limited domain sizes, of which only one is suitable as a template 
[TT] . Recently Xe buffer layer assisted self assembly of Co on a BN-layer forming a moire 
with Rh(lll) was achieved [H]. Though the placement of the Co clusters is regular, 
the method hampers from the very low temperature of fabrication, the limited thermal 
stability and the only incomplete filling of the template cells. As a final example, on 
the W(110)/C-R(15xl2) surface carbide arrangements of Au, Ag and Co clusters could 
be grown, the latter two with a high degree of order [151 IlH]- Though the thermal 
stability of the arrangements is rather good, no macroscopic extension of a superlattice 
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is achieved due to the persistent existence of terraces without the surface carbide. 

Graphene moires with noble metal surfaces are a new and unique support for cluster 
superlattices. This was shown first for graphene moires on Ir(lll) allowing Ir cluster 
superlattice formation [17]. Though a special material system, the binding mechanism 
is versatile. Through metal deposition on the graphene layer - which is only weakly 
interacting with the metal substrate in the absence of metals deposits [18] - graphene 
locally rehybridizes from sp^ to sp^ carbon bonds at a specific location in the moire unit 
cell and in between the substrate and deposit metal, thereby forming strong carbon metal 
bonds [19]. In between the substrate and deposit metal the sp^ rehybridized graphene 
has tetrahedral bond angles and may thus be considered as diamondlike. Within each 
moire unit cell rehybridization is possible where locally the carbon rings center around 
a threefold coordinated hep-site or a threefold coordinated fcc-site (both having three 
C-atoms on atop substrate sites). These locations are named hep region or fee region, 
respectively. Experimentally and by DFT calculations we find that Ir clusters bind by 
far stronger to hep regions [171 IlHl EQ] . 

Fig. la displays a cluster superlattice grown on a graphene flake formed by 
temperature programmed growth on Ir(lll) pi]. With this method graphene covers 
the surface only partially, which is suitable for scanning probe investigations, as it 
retains the metal surface partially for calibration of cluster size through island coverage 
in the graphene free areas. Also direct comparison of the properties of clusters and 
those of the deposit islands on the metal is possible. Recently, we employed chemical 
vapor deposition at high temperatures to grow graphene fully covering the substrate and 
displaying hardly any defects [211 122] • Using a dedicated combination of both methods 
[23] . we are not only able to ensure full coverage but also the graphene and Ir dense 
packed atomic rows to be parallel, i.e. no traces of rotational variants are present [24] . 

In this manuscript we will show that graphene moires as active templates are 
superior to other systems enabling the growth of two dimensional cluster superlattices 
on solid surfaces through a unique combination of properties: (i) the cluster binding 
mechanism is universal enabling growth of a large diversity of materials as superlattice; 
(ii) the superlattice order is extremely high with a completely filled lattice; (iii) if desired 
the superlattice extends macroscopically without uncovered substrate patches; (iv) 
fabrication is possible at room temperature; (v) the superlattices possess a reasonably 
high thermal stability and display absence of alloying and interdiffusion within a large 
temperature range; (vi) the cluster size is tunable and the size distribution is narrow. 

To demonstrate the universality of our approach in view of superlattice forming 
materials we primarily focussed on cluster materials with potentially interesting 
structural, magnetic, catalytic or optical properties. As heuristic guidelines for the 
suitability of a material to form a superlattice we considered three factors: (i) A large 
cohesive strength of the material as an indicator for the ability to form strong bonds, 
(ii) A large extension of a localized valence orbital of the deposit material allows it 
to efficiently interact with the graphene vr-bond and thus to initiate rehybridization 
to diamondlike carbon underneath the cluster, (iii) A certain match of the graphene 
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material 


strueture 


a/A 


rd/pm 


Ecoh/eV 


Ir 


fee 


2.715 


70.8 


6.94 


W 


bee 


2.741 


77.6 


8.90 


Re 


hep 


2.761 


73.9 


8.03 


Pt 


fee 


2.775 


65.9 


5.84 


Au 


fee 


2.884 


63.5 


3.81 


Fe 


bee 


2.483 


38.2 


4.28 


Ni 


fee 


2.492 


33.8 


4.44 



Table 1. Crystal structure, nearest neighbor distance a, valence d-orbital radius rd 
and cohesive energy i?coii of the tested materials for cluster superlattice formation on 
the graphene moire on Ir(lll). 

unit eell repeat distanee on Ir(lll) of 2.452A [20] and the nearest neighbour distance 
of the deposit material is neeessary to fit the first layer eluster atoms atop of every 
seeond C atom. As small clusters - which are the relevant sizes to start superlattice 
growth - have a smaller lattice parameter compared to bulk materials, and as Ir works 
perfectly as a eluster material ^7\, we considered 2. 7 A as an optimal nearest neighbor 
distanee [25]. Table compares the data for the tested materials. Tungsten was selected as 
a likely candidate for superlattice formation (all three figures are in favor of superlattice 
formation), but with a different crystal strueture than Ir. Re, Pt and Au were selected 
as as potentially interesting materials for catalysis. According to our guidelines we 
expected Re to be most likely a superlattice forming material, also Pt, but with a 
slightly smaller probability and Au as an unlikely candidate. Fe and Ni were selected as 
materials because of their ferromagnetism. According to their figure of merits, however, 
we did not expect superlattice formation. 

2. Experimental 

Experiments were carried out in an ultra high vacuum variable temperature scanning 
tunneling microscopy apparatus with a base pressure of 1 x 10~^°mbar equipped with 
a mass separated ion source, low energy electron diffraction, a gas inlet for ethylene 
exposure through a tube ending approximately 2 cm above the sample surface and a four 
pocket e-beam evaporator. The sample was cleaned by repeated cycles of sputtering with 
a beam of 2.5 keV Xe+ ions at 300 K and annealing to 1520 K. Graphene was prepared 
using temperature programmed growth (TPG) at 1470 K or chemical vapor deposition 
(CVD) at 1370 K [211 El]. TPG results in a graphene surface coverage of about 25% 
with flakes of a few 100 A to about 1000 A extension. Adding a CVD step results in 
up to full graphene surface coverage and domains of fim and larger extension. Cluster 
growth was performed through e-beam evaporation of well degased high purity metals 
with typical deposition rates of a few times 10~^ML/s where 1 ML (one monolayer) 
corresponds to the surface atomic density of Ir(lll). During deposition the pressure 
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Figure 1. STM topopgraphs of Ir(lll) with graphcnc flakes after deposition of 
an amount 6 of various metals at 300 K. (a) 6 = 0.20 ML Ir, average cluster size 
= 17atoms; (b) 6 = 0.25 ML Pt, Sav = 22 atoms; (c) 9 = 0.44ML W, 
Sav = 38 atoms; (d) 9 = 0.53 ML Re, Sav = 60 atoms; (e) 9 = 0.77 ML Fe, 
Sav = 420 atoms; (f) 9 = 0.25 ML Au, Sav = 100 atoms. Image size 700 A x 700 A. 

remained in the low 10~^°mbar range. Precise calibration of the deposited amount 
was obtained by STM image analysis of the fractional area of deposit islands on the 
clean Ir(lll) surface after a defined deposition time. If necessary the sample was gently 
annealed prior to coverage calibration. Through annealing the heteroepitaxial islands 
became compact and their size increased, thereby minimizing imaging errors associated 
with the finite STM tip size. Clusters were imaged at deposition temperature, if not 
indicated otherwise. However, for deposition above 300 K or after annealing the sample 
was quenched to 300 K for imaging. 

3. Results 

3.1. Cluster structure and superlattice formation at 300 K 

Figure 1 displays STM topographs after deposition of 0.2-0.8 ML on Ir(lll) partly 
covered with graphene flakes. The graphene flakes are typically attached to substrate 
steps, but also extend over one or several of them [Fig. 1(a), 1(b), 1(e)]. In the areas 
without graphene, deposit metal islands of monolayer height formed from the evaporated 
material. In Figs. l(c)-l(e) already second layer island nucleation took place. Depending 
on the deposited metal the island nucleation density on Ir(lll) varies considerably being 
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highest for the W and Re, the metals with the highest cohesive energy [Figs. 1(c) and 
(d)]. The deposit islands mostly reflect the threefold symmetry of the substrate. It 
is obvious from Fig. 1 that all deposited materials are pinned to graphene flakes to 
a certain extend and that graphene on Ir(lll) is in all cases much more sticky to the 
deposited metals than the surface of graphite [27] . However, not all materials form a 
cluster superlattice. 

Ir and Pt form superlattices of similar perfection [compare Figs. 1(a) and 1(b)]. 
For the represented O ~ 0.2 ML both materials exhibit two distinct height levels of the 
clusters indicating an out-of-plane texture of the cluster orientation. Distinct height 
levels are present also for larger and higher clusters up to the coalescence threshold 
(compare also Fig. 1 of reference |l7j). The apparent height differences between 3 ML 
and 4 ML clusters as well as those between 4 ML and 5 ML clusters are for Ir and Pt in 
the range of 2.2 A and 2.3 A, i.e. of the size of a monatomic step height hi on a (111) 
terrace [/ii,ir = 2.22 A and /ii,pt = 2.27 A]. The height differences in lower levels differ 
from these numbers. Specifically, the apparent height difference between the graphene 
substrate and the first cluster height level for Sav > 15 atoms is always found to be larger 
than 2.5 A while the difference between the first and second height level is typically below 
2.0 A. We interpret these deviations from the (111) step height as density of state effects. 
While for Ir and Pt the height of each single cluster is an integer number of (111) layers, 
the average cluster height /iav may be noninteger due to the averaging. Fig. 2 displays 
an analysis of h^v for Ir clusters grown at 350 K and Pt clusters grown at 300 K. The 
decrease of h^v for Ir clusters between 1.5 ML and 2.0 ML is due to the onset of cluster 
sintering. Upon cluster sintering the clusters reshape and material flows into the gaps 
separating the clusters, thereby causing frequently a height reduction. It is also apparent 
from Fig. 2, that large Pt clusters tend to grow flatter than Ir clusters, giving rise to a 
somewhat earlier cluster sintering. 

After having established the [111] out-of-plane texture of Pt and Ir clusters, the 
question arises whether the clusters possess also in-plane-texture, i.e. whether the 
clusters are also oriented within the surface plane. Fig. 3(a) displays a small area of an 
Ir cluster array after deposition of 1.5 ML imaged with a large tunneling resistance of 
2.3 X 10^ Q. Despite the absence of atomic resolution at least some cluster edges appear 
to be oriented along the (211)-direction. Lowering the tunneling resistance to 7 x 10^ Q 
- thereby bringing the tip close to the cluster surfaces - and tuning the contrast to the 
different levels of the cluster top mesas changes the picture. The top mesas are almost 
atomically resolved (some row like corrugation is visible) and the edges of the top mesas 
are unambiguously oriented along (110). It needs to be noted that under low tunneling 
resistance conditions necessary for atomic resolution clusters are usually picked up by 
the STM tip. The apparently different orientation of cluster edges at high tunneling 
resistances is a mere imaging artifact. If a tip scans a hexagonal grid of elevated objects 
in a large distance not the shape of the objects forming the grid but the hexagonal grid 
itself determines the apparent orientation of the boundaries between the objects. The 
inset of Fig. 3(d) represents a ball model of a five layer cluster containing 140 atoms 
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Figure 2. Ir (full squares) and Pt (red dots) average cluster heights hrmav in 
monolayers as a function of deposited amount 0. 




Figure 3. (a) STM topograph after deposition of 1.5 ML Ir on a grapliene moire on 
Ir(lll), Sav = 130. The dashed line indicates the position of an underlying substrate 
step. The visible clusters on the lower terrace have heights of 4 ML and 5 ML. For 
imaging the tunneling resistance was 2.3 x 10^ (/ = 0.5nAandJ7 = 1.15V. (b),(c),(d) 
Same location imaged successively with a low tunneling resistance of 7 x 10^ 
(/ = 30 nA and U = 0.2 V) and contrast tuned to the different cluster levels. Cluster 
edges are aligned along (110) substrate directions [see (b)]. Inset in (d): Ball model 
of a 5 layer cluster of 140 atoms consistent with the experimental observations. Image 
size 130 A x 160 A. 



which is consistent with the experiments. Although we are unable to prove directly 
that the cluster sidewalls are formed by {100} and {111} facets, facets of smaller slope 
would imply cluster contacts at their base. This is unlikely to be the case, as we find in 
annealing sequences clusters to reshape rapidly upon contact. To summarize, we have 
shown that Ir clusters (and the Pt ones most likely as well) are epitaxial clusters with 
the (111) cluster planes parallel to the substrate surface and the (110) cluster directions 
parallel to the (llO)-directions of the Ir substrate and the (li20)-directions of graphene. 
This epitaxy of Ir clusters on the graphene moire is also predicted by the geometry in 
the DFT based model of cluster binding [19] . 
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Also tungsten forms a cluster superlattice of high perfection. Compared to Pt and 
Ir the centers of mass of the clusters deviate slightly more from a perfect hexagonal 
superlattice. Still the scatter of the cluster positions is too small to result in cluster 
sintering for 6 = 0.44 ML as visualized in Fig. Ic. Close inspection of our data - also 
for = 0.04 ML - shows that less than 1% of clusters is out registry with the hep 
regions (see also below). The apparent cluster heights for = 0.44 ML range from 
4 A to 8 A with an average around 6 A. The clusters seem higher than the Ir ones for 
a comparable 0. Distinct height levels are also present for W clusters as is obvious 
from Fig. 1(c). However, the height levels are less well defined due to their smaller 
separation. Between adjacent clusters we measure frequently height differences of about 
0.8 A to 1 A. With a small distortion the dense packed W(llO) plane with a nearest 
neighbour distance of 2.74 A would fit to the triangular 2.46 A periodicity provided by 
the graphene for rehybridization. However, a [110] out-of-plane cluster texture would 
result in height levels spaced by 2.24 A, i.e. a similar spacing as for Ir and Pt. From 
the presence of intermediate levels we exclude this cluster orientation as the single one. 
Small step heights of 0.8 A to 1 A are consistent with a [111] out-of-plane texture with a 
step spacing of 0.91 A. The open W(lll) plane has also the threefold symmetry of the 
graphene moire in the hep or fee regions. However, in plane atomic spacings are 4.48 A. 
They could be fitted to the 4.27 A separation of a (a/3 x ^/3)R30 superstructure of 
graphene. For a final statement more detailed experiments are certainly necessary but 
also likely to be rewarding to settle the issue of the W cluster superlattice texture. 

For Re only partial order of the cluster arrangement is realized as visible in Fig. 
1(d). Two possible reasons could cause the apparent disorder: (i) the existence of 
several adsorption mimima within a moire unit cell, causing less ordered growth followed 
by early coalescence and cluster rearrangement; (ii) a too low depth of the cluster 
size dependent adsorption mimina. Assuming an increase of the potential well with 
increasing cluster size, adatoms and small clusters up to a certain size Sc would be 
likely to leave their unit cell during growth. Both scenarios would give rise to the 
observed heterogeneous size distribution, but the latter would explain also the existence 
of a large number of empty moire unit cells. For Re we performed also low coverage 
experiments [compare Fig. 7(d)]. We find a comparatively low moire unit cell occupation 
probability (or a comparatively large Sav) for the deposited amount, but with few 
exceptions all clusters adsorbed on a regular grid. We are therefore convinced that 
indeed the adsorption site minima are not of sufficient depth to keep small Re clusters 
up to a certain critical size (being much larger than the one for Ir, Pt and W) during and 
after growth on their position within the moire. We note that the poor quality of the 
Re cluster superlattice is at variance with our expectations, as Re has a large cohesive 
energy, an extended d-orbital and a reasonably well matching nearest neighbour distance 
(compare table 1). 

For Fe [Fig. 1(e)], Au [Fig. 1(f)] and Ni (data not shown) at room temperature 
no superlattice can be realized. Large clusters, lacking defined height levels or shape 
features are formed. The unstructured, hemispherical clusters display heights of up to 




Figure 4. STM topographs of Ir(lll) with graphcnc flalccs after depositing amounts Q 
of Re or Au at the indicated temperature T. (a) 6 = 0.45 ML Re, Sav = 41, T = 200 K; 
(b) O = 0.45 ML Re, Sav = 41, T = 200 K. Subsequently the sample has been annealed 
to 300 K and imaged; (c) 6 = 0.53 ML Re, Sav = 60, T = 300 K; (d) Occupation 
probability n of moire unit cells with clusters as a function of growth and annealing 
temperature T; triangles pointing to the right: Au clusters, 0.25 ML deposited at the 
indicated T; triangles pointing to the left: Au clusters, 0.25 ML deposited at 90 K 
and additionally annealed to the indicated T; down triangles: Re clusters, 0.45 ML 
deposited at the indicated T; up triangles: Re clusters, 0.45 ML deposited at 200 K 
and annealed to the indicated T. (e) 6 = 0.25 ML Au, Sav = 24, T = 90 K. 

23 A for Fe and 15 A for Au. The absence of a regular cluster superlattice for these 
materials is expected in view of their small cohesive energy and/or their limited valence 
orbital extension. Binding of the deposit metal to graphene is apparently too weak to 
trap adatoms and small clusters, i.e. the depth of the potential energy minima within a 
moire unit cell is not sufficient to stabilize a growing cluster. 

3.2. Low temperature cluster superlattice growth and annealing 

If we assume that cluster superlattice formation for Re and even more for Au, Fe and Ni 
is impeded by a too high mobility of small clusters during growth, lowering the growth 
temperature would decrease Sc and thus be an efficient strategy to improve superlattice 
formation. Fig. 4(a) displays a Re cluster superlattice grown at 200 K by deposition of 
G = 0.45 ML. Annealing to 300 K does not change the cluster superlattice as visible in 
Fig. 4(b). The side-by-side comparison to Re clusters grown at 300 K in Fig. 4(c) makes 
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the efficiency of lowering the growth temperature to improve superlattice formation 
obvious. The low growth temperature reduced s^, enabled the almost complete ffiling 
of the moire unit cells and allowed the clusters to grow to a size being stable even at 
300 K. One might speculate that a further lowering of the growth temperature would 
have improved the Re cluster superlattice even further. In the quantitative analysis of 
Fig. 4(d) it is also apparent that the difference in moire unit cell occupation between 
growth at 400 K and annealing to 400 K from a formation temperature of 200 K is even 
larger than the corresponding one at 300 K. At the same time it is also obvious from 
the data of Fig. 4(d) that growth at a low temperature and annealing preserves the 
superlattice only up to a limited temperature [(compare up and down triangles at 400 K 
in Fig. 4(d)]. Lowering the growth temperature to 200 K for Au did not significantly 
enhance cluster nucleation. However, lowering the growth temperature even further 
to 90 K resulted in an Au cluster superlattice of moderate order as shown in Fig. 4(e). 
Also in this case the once formed superlattice may be preserved to a higher temperature. 
However, as visible from the quantitative analysis shown in Fig. 4(d), the Au cluster 
superlattice is deteriorated already after annealing to 220 K. 

Mild annealing leads to more subtle effects which do not affect the positional order 
of the cluster array. For 0.45 ML of Ir deposited at 300 K, by annealing to 450 K for 
300 s the amount of single layered clusters decreases from (15 ± 2) % to (10 ± 2) % with 
a negligible decrease of the overall occupational density of the moire with clusters. This 
implies, that single layered clusters of a certain size transform to a more stable two 
layered form upon annealing. 

3.3. Cluster seeding 

Low cohesive energy metals tend not to form cluster superlattices on the graphene moire 
on Ir(lll), as tested by us for Au, Fe, and Ni. Such metals wet high cohesive energy 
metals, due to their lower surface free energy. High cohesive energy metals mostly form 
rather perfect cluster superlattices, tested here for Ir, Pt and W. It is thus natural 
to apply cluster seeding, i.e. to define the positions of the clusters by a small 9 of a 
high cohesive energy metal and to grow these seeds by subsequent deposition of a low 
cohesive energy metal [HI [12]. The successful application of this method is visualized 
in Fig. 5. As shown in Fig. 5(a) ffist through deposition of 0.1 ML Ir seed clusters 
are created in nearly all moire unit cells. Subsequent deposition of Au [Fig. 5(b)] or 
Fe [Fig. 5(c)] at 300 K results in highly perfect Au and Fe cluster superlattices with 
Ir cores. Comparing Fig. 1(f) with Fig. 5(c) for Au and of Fig. 1(e) with Fig. 5(d) 
for Fe makes the dramatic effect of seeding obvious. Even after deposition of amounts 
beyond the coalescence threshold the seeding has beneficial effects on cluster uniformity 
and distribution. For deposition of O = 2 ML Fe still a high density of large, uniformly 
sized and spaced Fe clusters grows. These Fe clusters contain 760 Fe atoms and have a 
height of 20 A [compare Fig. 5(h)]. They are not more positioned on a superlattice, but 
span 4 to 5 moire unit cells and contain about 40 Ir atoms. 
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Figure 5. (a), (c), (d),(c) STM topographs of Ir(lll) with graphcnc flakes after 
depositing amounts 8 of Ir and Au or Fe at 300 K. Image size 500 A x 500 A. 
(a) Q = 0.10 ML Ir. Sav = 9. (b) Line scan along the hne indicated in (a). The two red 
arrows indicate the height level of ML clusters, (c) <d = 0.10 ML Ir and subsequently 
= 0.35 ML Au corresponding to 30 Au atoms per seed cluster, (d) O = 0.10 ML 
Ir and subsequently Q = 0.70 ML Fe corresponding to 61 Fe atoms per seed cluster, 
(e) Q = 0.10 ML Ir and subsequently Q = 2.0 ML Fe. Average number of Fe atoms 
per cluster is 760. (f) Line scan along the line indicated in (d). (g) Line scan along 
the line indicated in (e). (h) Line scan along the line indicated in (f). The red dashed 
lines in (f)-(h) indicate the height level of monolayer Ir clusters in (b). 



3.4- Temperature stability 

For applications of cluster superlattices in nanomagnetism and nanocatalysis thermal 
stability of the cluster arrays and the absence of sintering at the temperature of use is of 
crucial importance. To provide data in this respect we investigated the thermal stability 
of the materials tested so far and display the results in Fig. 6. The example annealing 
sequence of Figs. 6(a) to 6(f) shows the gradual decay of the cluster superlattice through 
isochronal annealing steps of 300 s up to 650 K. The Pt cluster superlattice remains intact 
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Figure 7. (a)-(e) STM topographs after deposition of 0.01 ML Ir at 350 K and 
subsequent heating to 390 K, The time lapse sequence images the same surface spot 
every 120 s at 390 K. Image size 250 A x 250 A, Sav = 4.5 atoms. Circles indicate, where 
changes take place in successive images, (f)-(j) STM topographs after deposition of 
2.0 ML Ir at 350 K and subsequent heating to 470 K, The time lapse sequence images 
the same surface spot every 120 s at 470 K. Image size 150 A x 150 A. 

up to 400 K. Fig. M^g) quantifies annealing by plotting the temperature dependence of 
the moire unit cell occupation probability ?t, as a function of temperature T. It is 
apparent that the Ir cluster superlattice is indeed the most stable one, decaying as the 
only one in two steps. Most cluster superlattices are stable up to 400 K, which provides 
a reasonable temperature window for nanocatalysis and nanomagnetism experiments. 

The decay of all cluster superlattices occurs due to the thermally activated motion 
of clusters. The clusters fluctuate around their equilibrium positions within the moire 
unit cell. The magnitude of cluster fluctuations depends on the cluster size and the 
internal cluster structure (the isomer). Upon encounter during their fluctuations the 
clusters merge immediately, on a time scale of less than a second. The outcome of 
the merging again depends on cluster size. We distinguish two prototypical situations, 
(i) The clusters consist only out of a few atoms. Such clusters result from a very low 
deposited amount 6. For such a 6 also the moire unit cell occupation probability n is 
typically well below 1. Through thermally activated fluctuations a cluster may surmount 
the activation barrier Ea to leave its moire unit cell. If the cluster arrives in an empty 
cell it will rest there, if it arrives in an occupied cell the two clusters in the cell coalesce. 
They reshape completely such that as end product a single, compact cluster results 
which is located entirely within a single moire unit cell. This is the regime of complete 
cluster coalescence. Figs. 7(a)-(e) display a sequence of STM topographs taken at 390 K. 
The visible clusters result from O = 0.01 ML and have s^v = 4.5 atoms. White circles in 
subsequent stills indicate locations of thermally activated changes. The white circles in 
the upper left corner of Figs. 7(a)- (c) highlight a situation of thermally activated cluster 
motion resulting eventually in complete cluster coalescence. Note that the resulting 
cluster appears to be larger and higher in consequence of complete coalescence. In the 
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Pt annealing sequence of Fig. 6 in (d) a considerable number of larger and higher clusters 
appear, which are preferentially located next to empty moire cells. They are likely to be 
formed by complete coalescence. In addition, in Figs. 7(a)- (e) a number of thermally 
activated cluster jumps into empty cells are circled, which consequently do not result 
in cluster coalescence, (ii) If the clusters are large, close to the coalescence threshold, 
cluster merging proceeds differently. Still such clusters fluctuate around the location 
of their potential energy minimum. According to their large mass and better internal 
stability the magnitude of fluctuations has diminished and one might ask how cluster 
merging takes place at all. However, only small fluctuations are necessary to initiate 
cluster merging, as due to their large size the cluster have only a small edge separation 
of a few Angstroms. For coalescence the clusters do not have to leave their moire unit 
cell, but it is sufficient to move a little up their shallow potential energy depression to 
encounter a neighbouring cluster. Due to their size and significant binding the resulting 
cluster spans two moire unit cells and does not reshape completely. This is the regime 
of cluster sintering or incomplete cluster coalescence. Figs. [7] (f)-(j) display a situation 
of cluster sintering imaged at 470 K for clusters formed after deposition of 1.5 ML Ir 
(130 atoms per moire unit cell) at 350 K. Two subsequent cluster sintering events of 
neighbouring clusters result eventually in a new single cluster extending over three 
moire unit cells. The two scenarios depicted above are extreme cases and intermediate 
situations occur. From what has been said above it appears that cluster superlattice 
stability depends also on It is expected that arrays of medium sized clusters with 
diminished fiuctuation amplitudes and still sufficient separation from their neighbours 
are the most stable ones. 

To obtain a quantitative estimate for parameters determining cluster superlattice 
decay we model it as follows. We assume the cluster superlattice to consist of clusters 
with a unique activation energy for cluster interaction with a neighbouring cluster 
and an interaction frequency u = uqc'^'^^''^^ , where uq is an attempt frequency 
characterizing the frequency of cluster fiuctuations, the Boltzmann constant and T 
the temperature of the experiment. The probability that one cluster encounters another 
one is proportional to n. We assume complete coalescence, i.e. the final cluster to 
occupy only a single moire unit cell. The number of these events is as well proportional 
to n. Under these conditions the decrease of n with time t at a given temperature is: 



Using integration by parts we solve this differential equation for the boundary of the 
cluster density ni before, and n2 after an annealing step of a fixed time interval At 
resulting in 



We emphasize here that our crude approximation effectively averages over different 
size dependent interaction frequencies u for a given size distribution. However, using the 




(2) 
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65 Hz 


2700 Hz 


Ir, 0.45 ML 


0.28 eV 


0.08 eV 


0.06 Hz 


0.03 Hz 


0.05 Hz 


Pt, 0.25 ML 


0.60 eV 


0.08 eV 


500 Hz 


430 Hz 


3100 Hz 


Pt, 0.70 ML 


0.38 eV 


0.02 eV 


6.2 Hz 


2.3 Hz 


3.7Hz 


W, 0.44 ML 


0.47 eV 


0.04 eV 


33 Hz 


20 Hz 


52 Hz 



Table 2. Diffusion parameters and corresponding statistical errors as derived from 
v{T), for the cases of Ir, Pt an W. The parameters have as wcU been determined 
separately for two sections (I) and (II) of the dataset as indicated in Fig. [5] (h) Attempt 
frequencies have an statistical asymmetric error, so that negative and positive one are 
given. Systematic errors may be larger (see text). 

annealing time At and the annealing temperature T, this appoach allows one to derive 
the temperature dependence of v from an annealing sequence as shown in [6] (a)-(f). 

The resulting Arrhenius plots are shown in Fig. [7] (h). Activation energies are 
between 0.38 eV and 0.75 eV. The resulting pq lie between 1.4 Hz and 500 Hz. These 
attempt frequencies are much lower than a typical phonon frequency and also much 
lower than what is found for the diffusion of adatoms and small adclusters (compare 
e.g. [28]). According to transition state theory the low vq point to exceptionally large 
differences of the partition function of clusters in the bound state versus clusters in 
the transition state. The diffusion parameters have to be viewed as effective diffusion 
parameters for an ensemble of clusters comprising all size effects and not as the properties 
of an individual cluster. This is also illustrated by the large variations in parameters for 
the observed cluster lattices from Pt with different average cluster size. 

For the case of Ir as shown in Fig. [6] (g) and (h), there is a distinct discontinuity in 
the cluster density and consequently in the estimated interaction frequencies p between 
550 K and 650 K. Interestingly this transition coincides with all single layered clusters 
dying out. We interpret the discontinuity as a result of different diffusion parameters 
for single layered and multi layered clusters. 

This approach of effective diffusion parameters is checked for consistency with a 
kinetic Monte Carlo simulation. The algorithm for the simulation is based on work by 
Bortz at al. [26]. The cluster lattice is modeled as a hexagonal lattice with clusters, 
which can hop to adjacent sites with a frequency v = uqc^'^^^^'^ based on uq and 
as given in table [2l The cluster lattice dwells at each annealing temperature for 300 s 
and is heated successively to higher temperature annealing intervals. Time for cooling, 
imaging, and reheating is omitted. The simulation reproduces the cluster densities well 
as shown in figure [HI Kinks occur in the curve, where the temperature changes. Not 
surprising. Iridium is an exception to the good match because one kind of clusters dies 
out rapidly and the approximation of effective diffusion parameters breaks down for the 
examined temperature range. 
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Figure 8. Cluster density n evolution with a Monte Carlo simulation using average 
diffusion parameters. The temperature is increased every 300 s corresponding to Fig. 
Wis) ^-iid (h) . Cluster densities after the annealing steps as in figure [6] are rcshown for 
comparison. 

3. 5. Binding sites of clusters in the superlattice 

Experimentally [171 EHj and by calculations [19] we find Ir clusters to adsorb 
preferentially in hep regions of the moire unit cells. These sites differ from the fee 
regions only by the fact that instead of a threfold coordinated fee hollow site a threefold 
coordinated hep hollow site is centered in the carbon ring. They differ significantly from 
the atop-type area which has an atop site centered in the carbon ring. At low growth 
temperature we find for Ir deposition also fee regions to be populated by small clusters, 
consistent with the similarity of the two areas [20]. It is not at all evident that also other 
materials adsorb preferentially to hep regions. Specifically for materials with a different 
crystal structure like W with its bcc structure or Re with the hep crystal structure we 
would not be surprised to find the clusters adsorbed preferentially to fee regions or to 
be even unspecific to the small difference caused by the second layer underneath the Ir 
surface. To obtain cluster binding sites experimentally we make use of the fact that the 
graphene sheets on Ir form a jagged zig-zag edge when in contact with a (llO)/{100} 
microfacet or A-step of the substrate. The step undulation of the graphene sheet has 
the moire periodicity. The protrusions of the graphene flake's edge are bowing out 
towards the Ir terrace at atop-type areas. This fact does evidently not depend on the 
deposited material and allows us unambiguous cluster adsorption site assignment. In 
Fig. [9] for Ir, Pt, W and Re the corners of the moire unit cell grids are fixed to these 
atop-type areas. We find that the clusters are always located in the triangular half-unit 
cells pointing away from the one dimensional graphene-Ir interface (green triangles in 
Fig. [9]). According to our unit cell assignment (compare Fig. 1 of reference [17]) these 
are hep regions. 




Figure 9. Binding site determination for (a) Ir. (b) Pt. (c) W and (d) Re clusters. 
Deposition was performed at 300 K and (a) Q ~ 0.10 ML, Sav = 9 atoms; (b) 
e = 0.04 ML, Sav 6 atoms; (c) 6 = 0.04 ML, Sav = 4 atoms; (d) 6 = 0.03 ML, 
Sav — 10 atoms. Due to the moire the graphene flakes form a jagged edge if in contact 
with a (110)/{100} microfacet or A-step. The dark tips of the jagged graphene edge 
are atop site areas. Fixing the grid of moire unit cells to these positions enables a 
binding site assignment. The clusters always sit in the green triangular half- unit cells 
pointing away from the step (see text). 

Let us point out that in Fig. 8 of reference [20] and in the accompanying text of 
Sec. 8 the labehng and use of hep region and fee region is erroneously interchanged. 

3. 6. Towards cluster superlattice materials 

To probe the properties of cluster superlattices by averaging techniques and to 
investigate their suitability for potential applications it is necessary to cover a sample 
macroscopically with a cluster superlattice. This need is evident if one considers e.g. the 
analysis of reaction products from a cluster superlattice in nanocatalysis. The presence 
of the bare metal would result in additional peaks in thermal desorption spectra and 
certainly complicate the data analysis. To establish a macroscopic cluster superlattice 
the graphene moire must cover the metal substrate entirely, be of unique orientation 
with an as large as possible moire supercell which displays upon deposition of suitable 
materials local rehybridization. While it is likely that a number of graphene moires 
on different metals fulfill all these conditions, so far they have been proven only for 
graphene moires on Ir(lll). Some optimization of the graphene growth procedure was 
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Figure 10. Ir cluster superlattice grown at 300 K with 8 = 0.80 ML resulting 
in Sav = 70 on graphcnc prepared by temperature programmed growth followed 
by chemical vapor deposition and extending over the entire sample. Image size 
0.5 ^m X 0.3 ^m, inset 500 A x 300 A. 

necessary to achieve simultaneously full graphene coverage and a single orientation of 
the graphene and the graphene moire [21], [221 123] • FiS- 113 displays a large scale STM 
topograph visualizing to a certain extent the quality of the available substrate. The 
cluster superlattice is present in the entire topograph in unique orientation and even 
steps merely present locations where a line of clusters is missing, but without disturbing 
the overall alignment of the superlattice. 

4. Conclusion 

In conclusion, we have established that the graphene moire on Ir(lll) is a versatile and 
active template for cluster superlattice growth of a great variety of materials and with 
macroscopic lateral extension. If necessary, techniques like low temperature growth or 
cluster seeding may permit cluster superlattice growth for cases, where simple room 
temperature deposition fails. The high thermal stability of the cluster superlattices and 
the ability to grow them on macroscopic areas opens new opportunities for fundamental 
cluster research and applications. 
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